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ABSTRACT 

Using stellar dynamics and strong gravitational lensing as complementary probes. ISand et alJ {2002, 
2003?) have recently claimed strong evidence for shallow dark matter density profiles in several lensing 
clusters, which may conflict with predictions of the Cold Dark Matter paradigm. However, systematic 
uncertainties in the analysi s weaken th e constraints. By re-analyzing their data, we argue that the 
tight constraints claimed bv ISand et aLl were driven by prior assumptions. Relaxing the a ssumptions, 
we find that no strong constraints may be derived on the dark matter inner profile from the lSand et alJ 
data; we find satisfactory fits (with reasonable parameters) for a wide range of inner slopes p oc 
with < < 1.4. Useful constraints on the mass distributions of lensing clusters can still be obtained, 
but they require moving beyond mere measurements of lensing critical radii into the realm of detailed 
lens modeling. 

Subject headings: gravitational lensing — dark matter — galaxies: clusters — galaxies: clusters: 
individual (Abell 383, MS 2137-23) 



1. COLD DARK MATTER OR CORED DARK MATTER? 

The central density profile of dark matter halos has 
been the subject of considerable debate in recent years. 
On the theoretical side, numerous groups have claimed 
that dissipationless N-body simulations in the Cold Dark 
Matter (CDM) model produce universal halo profiles. 
Despite vigorous debate over the exact shape of the in- 
ner profile, there is general agreement that predicted 
CDM halos hav e central density cusps, p oc with 
~ 1-1.5 (e.g-. lMoore et alJll99Sl iNavarro et al.lll997l 
2003). On the observational side, similar controversy 
has raged over the question of whether such cusps are 
present in real galaxies, and whether their absence would 
challenge the CDM paradigm (e.g., iDutton et al]l2003T: 
iSimon et alJ ^OOl^. and references therein). The most 
difficult issue is accounting for the effects of baryonic 
matter, which contributes directly to the gravitational 
potential an d may also modify the dark matter distri- 
bution (e.g..lBlumenthal et alJll986tlEl-Zant et alJl2003t 
lLoeb fc Peeblesll2003|) ~ 

Strongly lensed arcs in galaxy clusters probe the grav- 
itational potential on scales (r ~ 50-100 kpc) large 
enough to avoid significant baryonic contamination, and 
hence can provide a relatively clean test of the predictions 
of dissipationless N-body simulations. Several groups 
have constrained the density profiles of i ndividual clus- 
ters b y studying their lensing properties lA.th.reva et afl 
120021 IClowe Zr. Schneider! 12002 IDahlel 120031) . and Tn 
some cases have obtained stringent limits from combined 
stron g- and weak -lensing analyses (e.g., IGavazzi et al.l 
l2^0llKneib et al.ll2003^l . 

Another approach is to disentangle the baryonic and 
dark matter contributions to the net potential on small 
scales, similar to studies of galaxy rotation curves. 
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IKelson et alJ l|2002j) have used the velocity dispersion 
profile of the cD galaxy in the cluster Abell 2199 to- 
gether with the kinematics of the cluster members to de- 
compose the stellar and dark matter components. They 
found that a = 1.5 dark matter cusp is ruled out by the 
velocity dispersion data, while a = 1 cusp is difficult to 
reconci le with a reasonab l e mas s-to-light ratio for the cD 
galaxy. ISand et all l)2002l 120031 hereafter S03) have gone 
a step further and combined dynamical and strong lens- 
ing analyses for six clusters. Two of their clusters provide 
particularly tight and interesting results: inner slopes of 
= 0.57 ±0.11 for MS 2137-23 and = 0.38 ±0.06 for 
Abell 383, which strongly conflict with the steep cusps 
predicted for CDM. Such a conflict on cluster scales may 
be more troubling than similar conflicts on galaxy scales, 
as there are fewer proposed processes that are capable of 
disrupting dark matter cusps in clusters. 

The S03 results are somewhat puzzling, however, as 
illustrated in Figure ^ This figure shows velocity dis- 
persion profiles and representative lensed arcs for three 
models of Abell 383: the best-fit model found by S03, a 
model with a = 1 dark matter cusp, and a model with 
= 0. While the fit with a = 1 cusp is by no means 
perfect, it is clearly not ruled out at the 10-cr level (as 
claimed by S03). Moreover, it is not significantly worse 
than the best-fit model presented by S03. 

In this paper we study the origin of this discrepancy. 
We argue that the tight constraints claimed by S03 re- 
sult from assumptions made in their modeling. In par- 
ticular, they assumed spherical symmetry and adopted 
particular values for the scale radius at which the dark 
matter density profile transitions from r~ l3 to r~ 3 . Re- 
laxing these assumptions, we find that no interesting con- 
straints on the dark matter profile can be derived from 
the S03 data alone. However, more detailed modeling of 
the lensing properties of these clusters may yield stronger 
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Fig. 1. — (a— c) Representative model arcs for Abell 383. The 
model parameters are listed in Table [I] The solid curves show 
the arcs, while the dotted curves show the lensing critical curves. 
Model a is the best-fit model from S03; the second tangential arc 
is an artifact of the assumed spherical symmetry. In models b and 
c, the arclet near the tangential arc is a predicted counter-image of 
the radial arc. (d) Velocity dispersion profiles for models a (solid 
line), b (dashed line), and c (dash-dot line), compared with the 
data from S03. 



Table 1. Sample Model Parameters 



Panel 


a 


b 


c 


M 9 [M ] 


1.8 x 10 12 


10 12 


1.2 x 10 12 


na 


13.75" 


13.75" 


13.75" 


eg 





0.2 


0.2 


P S 


2 


2.2 


2.2 




1.3 x icr 2 


3.8 x 10" 3 


4.2 x 10~ 3 


M h [M Q ] 


4.6 x 10 15 


1.9 x 10 15 


4.9 x 10 14 


C v j r 


9.6 


4.7 


24 


rh 


127" 


194" 


24" 


<'h 





0.2 


0.15 


Ph 


0.4 


1 







0.576 


0.166 


1.33 



Note. — Parameters for the models of Abell 383 shown in Fig- 
ure [I] We use zi = 0.189, z a = 1.0 in a flat ACDM cosmology 
with Qm = 0.3 and h = 0.7. Quoted halo masses are virial masses, 
where t he virial overdensity is ap proximated using the fitting for- 
mula of Bryan & Norman (1998). Concentrations are defined by 
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constraints. 

After the completio n of this paper, we became aware o f 
independent work bv iBartelmann fc Meneghettil |2003), 
who reach similar conclusions. 

2. ANALYSIS 

2.1. Data and models 

The observables modeled by S03 are the lensing critical 
radii (basically the radii of the tangential and radial arcs) 
and the velocity dispersion profiles of the brightest clus- 
ter galaxies; the data are given in their Tables 4-5. For 



a given mass model, the critical radii are determined by 
projecting the density distribution and locating the sin- 
gular points of the lens mapping (see e.g. lSchneider et all 
1992), while the velocity dispersion profile is computed 
by solving the spherical Jeans equation (see appendix) . 

We use a generalization of the two-component mass 
model used by S03. For the stella r component, we use 
the r] model llTremaine et al.lll994|) . 
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S03 mainly used a Jaffe model, which is a particular case 
of the i] model with (3 g — 2. Following S03, we describe 
the dark matter with a generalized NFW-type profile, 
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2.2. 



A case study: Abell 383 

In the S03 results Abell 383 provided the strongest 
evidence that the dark matter slope is shallower than 
so we use this cluster to study important system- 
atic uncertainties in the analysis. If we center both the 
stellar and dark matter models on the observed central 
galaxy position, orient them along the observed galaxy 
position angle, and fix the stellar model's half-light radius 
to the observed galaxy effective radius (13.75"), then the 
remaining model parameters are: the galaxy mass, in- 
ner slope, and ellipticity (M g , /3 g , e g ); and the halo mass, 
scale radius, inner slope, and ellipticity (Mh,rh, Phi^h)- 
S03 used a spherical Jaffe model for the stellar compo- 
nent (f3 g — 2 and e g = 0), and a spherical halo (eh = 0) 
with scale radius = 400 kpc. With the same assump- 
tions, we reproduce their constraints on the dark matter 
inner slope (3h- 

The sur face brightn e ss pr ofile of the central galaxy, 
shown by Smi th et al.l l|200H) . is steeper than expected 
for a Jaffe model and so we adopt (3 g s» 2.2, however 
we note that this steep slope has been disputed (R. Ellis 
2003, private communication). We find that with this 
stellar component the models can accommodate a some- 
what steeper dark matter profile: with a f3 g = 2 galaxy 
the best-fit halo has (3h ~ 0.38, while with j3 g — 2.2 it 
shifts to flh ~ 0.45. The shift of A(3h 0.1 is compara- 
ble to various systematic effects considered by S03. The 
model with a g = 2.2 galaxy fits better than the S03 
model (with A\ 2 — —1-5), so henceforth we focus on it. 

A larger effect is associated with the dark matter scale 
radius. The S03 constraints on (3h depend crucially on 
their assumption of = 400 kpc, because there is a de- 
generacy between rh and the dark matter slope (3h- For 
example, assuming rh = 200 kpc would yield j3h ~ 0.18, 
while assuming rh = 800 kpc would yield (3h ~ 0.66 
(for spherical models). The latter model fits considerably 
better than the S03 model (Ax 2 = —2.6), even though 
S03 claimed that 0h > 0.55 was excluded at 99% confi- 
dence. For spherical models, larger [3h do require large 
scale radii r^: for example, the best-fitting model with 
0h = 0.8 has rh > 1 Mpc! Such implausibly large scale 
radii (or small concentrations) would suggest that steep 
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Fig. 2. — Constraints on dark matter inner slope f)^ and ellip- 
ticity eh for Abell 383. Contours are drawn at the 68, 90, 95 and 
99% confidence levels. 



cusps are still disfavored, but that is very different from 
claiming that steep cusps are inconsistent with the data. 
The lesson is that must be treated as a free parameter, 
while keeping in mind that its value should be checked a 
posteriori for plausibility. 

The most dramatic effects arise upon dropping the as- 
sumption of spherical symmetry. The ellipticity of the 
galaxy is not exceedingly important; we use the observed 
value e g = 0.2 and find that it does not significantly 
change the results. What matters most is the elliptic- 
ity of the halo. Figure [21 shows likelihood contours in the 
plane of eh and Ph, optimizing over the other parameters 
(M g , Mh, and rt). We find that a broad range of dark 
matter slopes are consistent with the data. In the limit 
of spherical symmetry we recover the tight constraints 
on j3h found by S03, but if we allow even a relatively 
small ellipticity we find successful models over the range 
< < 1.4. Thus, the limits found by S03 appear 
to be an artifact of their prior assumptions for the mass 
model, notably the assumption of spherical symmetry. 
Note that the models we find with steeper inner slopes 
(/3 rj 1) have perfectly sensible parameters from a the- 
oretical standpoint. For example, model b in Table ^ 
with a Ph = 1 cusp, has galaxy mass M g = 10 12 Af Q , 
halo mass Mh = 1.9 x 10 15 Mq, and a halo scale ra- 
dius Th = 610 kpc, corresponding to a concentration of 
c V ir = 4.7. These halo parameters are fully consistent 
with the virial masses and concentrations expected for 
hot (T = 7.1 keV) X-ray luminous clusters. (Note that 
with moderate ellipticities it is possible to have Ph = 1 
without an unphysically large scale radius.) 

2.3. Additional constraints: MS 2137-23 

The case of Abell 383 suggests that simply combin- 
ing the lensing critical radii with dynamical data cannot 
strongly constrain the dark matter slope. Fortunately, 
the use of more detailed lensing data and modeling can 
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Fig. 3. — Likelihood as a function of the dark matter slope (3^ in 
MS 2f37— 23, based on lensing, X-ray temperature, and dynamics. 



provide constraints on the ellipticity that significantly 
improve the constraints on the dark matter slope. To 
illustrate, we consider MS 2 137-23. While S03 used 
only the lensing critical radii. IGavazzi etaD p003l) iden- 
tified multiple images of 26 distinct sources in the arcs 
produced by this cluster, enabling much more detailed 
modeling. Additionally, the X- ray temperature h as been 
measured to be T = 5.56 keV {Allen et alJl2001|) . which 
can provide constrai nts on the halo mas s via the M-T re- 
lation determined bv lAllen et alJ <|2001|) for relaxed lens- 
ing clusters. We combine these data with the velocity dis- 
persions from S03, but we inflate the positi onal error bars 
on the multiply imaged knots reported bv IGavazzi et alJ 
to 1", and inflate the error on the X-ray temperature to 
1 keV, in order to maximize the impact of the velocity 
dis persions on the fit. (We neglect the 5th image claimed 
bv IGavazzi et all as its detection is tentative.) 

We again use two-component mass models, assuming 
the galaxy and dark matter halo to be concentric. We 
use the observed values of the galaxy half-light radius 
(5.02"), ellipticity (e g = 0.17), and position angle. If we 
force the halo to be spherical and to have scale radius 
r/j = 400 kpc we recover the same constraints on Ph as 
S03. However, if we allow the dark matter halo param- 
eters to vary and optimize over them (and also over the 
galaxy mass M g ), we find the likelihood as a function of 
the dark matter slope Ph shown in Figure 03 The dark 
matter slope is constrained to be Ph = 1±0.35 (95% con- 
fid ence). This is quite consistent with the results found 
by IGavazzi et all which perhaps is not too surprising: 
since we found earlier that the velocity dispersion data 
are not terribly restrictive, it is reasonable that the joint 
constraints from dynamics and detailed lens modeling 
are similar to those obtained from lens modeling alone. 
Typical halo virial masses and concentrations obtained 
were roughly Mh ~ 7 x 10 14 Af Q and c v „ ~ 7. Inciden- 
tally, we find that mild halo ellipticities (eh ~ 0.2) and 
mild misalignment between halo and galaxy (AO < 10°) 
are favored by the fits. 

3. DISCUSSION 

We have argued th at the stringent constraints claimed 
bv lSand et alJ l)2003j) on the inner slope of the dark mat- 
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ter profiles in clusters — @h — 0.38 ± 0.06 for Abell 383, 
f3 h = 0.57±0.11 for MS 2137-23 — are not supported by 
their data. We have found successful models with dark 
matter slopes well outside these bounds. The models 
with steep slopes (f3h ~ 1) have sensible parameters: halo 
virial masses Mh ~ 10 15 M Q and concentrations c v ; r w 4- 
5. 

It appears that the tight constraints obtained by S03 
on the inner dark matter profile are artifacts of several 
simplifying assumptions in their modeling analysis. We 
have shown that the assumption of spherical symmetry 
and of a particular value of the halo scale radius artifi- 
cially restricts the range of parameters consistent with 
the data. It may seem surprising that small departures 
from sphericity (e.g., ellipticities e ~ 0.1-0.2) generate 
such striking differences in the lensi ng properties of the 
clus ters. However, as discussed bv iDalal et al.l (l2003f) 
and iBartelmann fc Meneghettil l)2003|) . the strong lens- 
ing cross section for objects with shallow radial profiles is 
highly sensitive to small shear perturbations, arising for 
example from ellipticity. In other words, small changes 
in ellipticity elicit large changes in the size of the critical 
curves for lenses with shallow profiles. This extra degree 
of freedom in the models thwarts any attempts to derive 
tight constraints on the density profiles from the critical 
radii alone. 

The lack of strong constraints does not mean that 
lensed arcs are not useful probes of the inner structure 
of cluster dark matter halos. We have also shown that 
detailed modeling of lensed features and dynamics does 
allow interesting constraints to be placed on the inner 
potential. Decomposing the central mass distribution in 
MS 2137—23 into stellar and dark matter components, 
we found constraints on the dark matter inner slope of 
= 1 ± 0.35 at 95% confidence. 

We wouldn't take this measurement too seriously, 
though. First, we have not fully explored the model pa- 
rameter space — considering tidal shear from nearby ha- 



los or substructures within the cluster, or a radially vary- 
ing ellipticity, could weaken the bounds on (3h- Second, 
it can be argued that decomposing the total mass into 
stellar and dark matter components, which is a standard 
approach, is not ideal. Because the effects of baryonic 
mass on the dark matter distribution are poorly under- 
stood, the dark matter profile obtained by subtracting 
the baryonic component from the total mass might be 
very different from the profile the dark matter would have 
assume d in the absence of bary onic interference. For ex- 
ample, lLoeb fc Peebles! l)2003|) have suggested that the 
decomposition procedure could give the appearance of 
a cored dark matter distribution, even for proper CDM 
halos. 

Thus, as with galaxy rotation curves, probing the dark 
matter on scales dominated by baryons is a perilous ex- 
ercise. A much better test of CDM predictions would be 
to measure the profile on scales that are DM dominated. 
Fortunately, giant arcs offer such a probe, since they oc- 
cur on scales (~ 50-100 kpc) beyond the influence of cen- 
tral galaxies, yet well inside the p oc regime expected 
in clusters. Simple measurements of lensing critical radii 
are not enough, however; detailed lens modeling is re- 
quired. Fortunately, the data to support such modeling 
are available, and several clusters have already yielded in- 
teresting (and in some cases surprisi ng) measurements of 
their dark matter distributions (e.g., lC.Mvnz/i et ^MM. 
IKneib et al]l200^ . 
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APPENDIX 

Here we discuss the calculation of the velocity dispersion and lensing properties of the mass models defined in §2.1. 
For spheri cally symmetric distributi ons with isotropic velocity dispersion tensors, the spherical Jeans equations are 
solved by ( Bi nnev fc Tremainel H987'l 

a V) = ~~r~\ / P( r ) — 7T^ dr > (!) 



-dr', 

p(r) J r ry ' r' 2 

where p(r) is the stellar density, G is Newton's constant and M(r) is the total (stellar+DM) mass interior to radius 
r. This velocity dispersion is not directly observed; instead the projected, luminosity-weighted dispersion is measured. 
Assuming a constant mass-to-light ratio, this becomes 

$™ P{r)a 2 (r) ^ r J_ RI dr 2 
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S(-R) J R 



R 2 dr 



where S(r) is the projected surface density. Averaging over a finite radial bin from R\ to R2 gives 

2 f „ j£o>(WWdR hf£ftr)GM<r)[Fft)-F&)]dr 

a hin {Hi,H 2 ) - 



/£ Z(R)dR f£ rp(r) [A (&) - A (Sl)] dr 

where 

A(x) = { l<\ and F(x 

For the density profiles we have employed, of the form 

p(x = r/r s ) = 



(2) 



(3) 



aVl - x 2 + sin 1 (x) - | x < 1 
x > 1 
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(4) 
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there are no closed-form expressi ons for the projected surface density in terms of elementary functions. The convergence 
K = S/Scrit takes the form (e.g., IWvithe et alJl2001fl 

«(« = R/r.) = 2k.u x ^ C' 2 {S ™ 6)n ~le dO, (5) 
Jo (u + smO) n p 

where k s — p s r s /T, cr i t and as usual X crrt is the lensing critical surface density. For elliptical surface density profiles 
with axis ratio q = 1 — e, we use k(^J x 2 + y 2 /q 2 ) with k s — p s r s / (gE cr it)- iKeetonl (|20D If) discusses how the deflection 
angle and distortion tensor may be expressed as one-dimensional integrals over k and its derivatives. The ellipticities 
of interest to us are small enough (e ~ 0.1-0.2 in the density, even smaller in the potential) that we assume it 
is a reasonable a pproximation to use the spherical Jeans equations to compute the line-of-sight velocity dispersion 
l|Kochaneklll99l . 
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